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Introduction 

Don  W.  Kinsey 
CFD  Research  Branch 


This  technical  report  was  produced  by  the  Computa¬ 
tional  Fluid  Dynamics  Branch  of  the  Aeromechanics 
Di'  ision,  Flight  Dynamics  Directorate,  Wright  Labora¬ 
tory  The  CFD  Branch  is  composed  of  two  Sections;  the 
CFD  Research  Section,  and  the  Interdisciplinary  and 
Applied  CFD  Section  The  CFD  Research  Section  con¬ 
ducts  basic  research  which  is  largely  sponsored  by  the 
Air  Force  Office  of  Scientific  Research.  The  lACFD  Sec¬ 
tion  conducts  research  and  development  of  a  more  ap¬ 
plied  nature.  Thus  the  CFD  Branch  as  a  whole  covers 
a  very  wide  range  of  technologies  from  geometry  gener¬ 
ation  to  post-processing  of  CFD  results. 

CFD  technology  is  expanding  very  rapidly.  New  algo¬ 
rithms,  techniques  and  applications  are  being  added  al¬ 
most  daily.  Some  of  these  represent  significant  advances 
in  the  state-of-the-art,  while  others  are  of  value  to  a 
more  limited  audience.  Most  of  the  more  complete  re¬ 
sults  are  published  in  symposium  papers  and  journal 
articles,  but  many  smaller  efforts  are  never  published. 
In  either  case  it  is  easy  for  the  work  to  become  lost:  the 
published  work  is  lost  in  the  volumes  of  papers  and  ar¬ 
ticles  that  we  never  seem  to  have  time  to  read,  and  the 
unpublished  works  are  lost  in  the  files  of  the  researcher. 
The  objective  of  this  report  is  to  document  in  a  single 
source  the  preceding  year's  CFD  accomplishments. 


The  technical  descriptions  provided  m  this  report  fol¬ 
low  closely  the  format  used  in  N.AS.A  s  National  .Aero¬ 
dynamic  Simulation  annual  report.  The  articles  are  de¬ 
signed  to  provide  enough  information  to  clearly  define 
the  objective,  approach  and  significant  accomplishments 
of  the  task,  but  be  short  enough  to  hold  the  reader  s  in¬ 
terest.  Where  applicable,  a  reference  that  contains  more 
detailed  information  is  provided.  In  all  cases,  the  name 
of  the  author  is  provided  and  he  may  be  contacted  di¬ 
rectly  for  details  about  his  work.  .An  alphabetic  listing 
of  each  author,  his  phone  number,  mailing  address  and 
e-mail  address  is  provided  on  the  last  of  the  report. 

The  picture  chosen  to  go  with  this  introduction  is  rep¬ 
resentative  of  our  current  capability  in  applied  CFD. 
The  geometr}'  was  digitized  from  blueprints,  the  dif¬ 
ferent  parts  of  the  geometry  were  combined,  and  the 
surface  grid  and  block  faces  were  produced  using  our 
I.3G/VIRGO  interactive  graphics  package.  The  volume 
(3-D)  grid  was  computed  with  a  code  called  PLUTO.  .An 
Euler  solver  called  MERCURY,  based  on  Jameson's  fi¬ 
nite  volume  technique,  produced  the  flow  solution.  The 
surface  geometry  and  particle  traces  were  produced  with 
the  .NASA-developed  PLOT3D  graphics  package. 


Euler  Solution  on  a  Complex  Fighter  Geometry 
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Three-Dimensional  Unstructured  Grid  Generation 

Frank  C.  Witzeman  Jr 
Interdisciplinary  and  Applied  CFD  Section 


Research  Objective 

One  of  the  most  difficult  and  time-consuming  tasks  in 
the  CFD  solution  process  is  grid  generation.  Interac¬ 
tive  graphics- based  programs  such  as  GRIDGEN  and 
I3G/VIRGO  have  overcome  many  of  the  limitations  and 
labor-intensive  tasks  of  multiple-block  structured  grid 
generation  .\lthough  these  methods  (and  many  others) 
have  enabled  CFD  solution  turnaround  times  on  the  or¬ 
der  of  weeks  or  days,  they  cannot  meet  design  environ¬ 
ment  goals  of  days  or  hours.  This  effort  is  focused  on  re¬ 
ducing  CFD  solution  times  from  weeks  to  days  through 
the  development  of  a  highly-automated,  easy-to-use  un¬ 
structured  grid  generation  procedure.  This  method  will 
be  applicable  to  arbitrary  and  complex  vehicles  and 
components  vvhich  are  nearly  impossible  to  model  with 
the  structured  grid  methods. 

Approach 

The  unstructured  grid  generation  method  is  part  of  the 
overall  Multi-Bodied  CFD  Method  being  developed  un¬ 
der  contract  to  the  Computational  Mechanics  Company 
(Dr  S.  R.  Kennon,  Principal  Investigator).  The  method 
is  based  on  a  generalized  three-dimensional  Delaunay 
concept  which  produces  grids  comprised  of  tetrahedral 
elements  instead  of  structured,  hexahedral  elements. 
For  most  cases,  the  only  input  required  is  an  accurate 
description  of  the  geometry.  Currently,  this  geometry 
description  is  the  same  as  that  used  for  the  structured- 
grid  methods.  Unique  features  of  this  method  are  the 
automatic  addition  of  grid  points  (both  on-  and  ofF- 


body)  and  the  enforcement  of  triangular  faces  which  lie 
on  the  geometric  surfaces.  Further  advantages  of  the 
method  include  the  automatic  insertion  of  grid  points 
in  regions  of  interest  (adaptive  gridding)  and  relocation 
of  grid  points  for  moving-body  simulations 

Accomplishments 

A  preliminary  version  of  the  unstructured  grid  genera¬ 
tor.  TETMESH.  has  been  produced.  The  software  is  be¬ 
ing  evaluated  by  several  Air  Force  organizations.  Three- 
dimensional  unstructured  grids  have  been  constructed 
for  several  test  cases,  and  Euler  solutions  have  been  ob¬ 
tained  using  our  in-house  flow  solver  COBALT,  For  the 
example  illustrated  below,  approximately  1  week  was 
required  to  generate  a  multiple-block  structured  grid. 
For  the  same  configuration,  only  6  hours  were  needed  to 
generate  the  unstructured  grid. 

Future  Plans 

The  interactive  three-dimensional  unstructured  grid 
generator  vvill  be  modified  to  be  more  robust  and  user- 
friendly.  More  generalized  criteria  will  also  be  developed 
for  unstructured  grids  made  up  of  irregular  tetrahedral 
and  triangular-prismatic  elements.  These  types  of  grids 
are  required  for  accurate  simulation  of  viscous  flows. 

Publications 

Kennon.  S.  R..  et  al.  'Geometry  Based  Delaunay  Tetra- 
hedralization  and  Mesh  Movement  Strategies  for  Multi- 
Body  CFD."  AIAA  Paper  9'2--t57.5.  August  1992. 


Navier-Stokes  Simulation  on  Unstructured  Hexahedral  Meshes 

Captain  Michael  J.  Aftosmis 
CFD  Research  Section 


Research  Objective 

This  research  effort  develops  and  evaluates  an  adaptive, 
high-resolution.  upwind  algorithm  for  the  simulation  of 
iiiviscid  and  viscid  flows.  Upw'ind  methods  provide  dis¬ 
crete  representations  of  shocks  and  contact  discontinu¬ 
ities  using  as  few  cells  as  possible,  and  adaptive  mesh 
techniques  automatically  adjust  the  local  flow  physics. 
The  combination  of  these  two  approaches  suggests  a 
natural  and  extremely  flexible  platform  for  simulating 
flowfields  with  features  spanning  a  variety  of  disparate- 
length  scales.  This  work  seeks  first  to  overcome  .some  of 
the  practical  and  implementational  issues  involved  with 
the  development  of  such  a  technique,  and  then  to  evalu¬ 
ate  the  utility  of  the  resulting  method  in  its  application 
to  a  variety  of  problems  of  interest. 

Appi’oach 

This  effort  applies  a  Roe-based,  finite-volume,  upwdnd 
technique  on  an  adaptive  mesh  consisting  of  unstruc¬ 
tured  hexahedral  cells.  Viscous  terms  are  modeled  with 
compact  central  differences.  As  the  solution  evolves,  the 
mesh  cells  adapt  through  directional  ceil  division  in  or¬ 
der  to  increase  the  resolution  of  the  computation  in  and 
around  emerging  flow  features.  By  concentrating  new 
computational  nodes  only  where  the  solution  requires 
further  enhancement,  the  method  uses  far  fewer  cells 
than  comparable  structured  mesh  techniques. 

Accomplishments 

An  unstructured  computer  code  was  written  and  vali 


dated  on  a  variety  of  inviscid  and  viscous  test  cases.  Nu¬ 
merical  simulations  of  3D  supersonic  corner  flows,  sub- 
aiid  transonic  airfoil  flows  and  simulation  of  laminar 
flow  over  a  delta  wing  were  all  in  ^reement  with  exper¬ 
imental  and  theoretical  results  The  method  typically 
requires  10  to  20  times  fewer  cells  than  an  e(|uivaiently 
resolved  solution  on  a  structured  mesh.  Moreover,  since 
the  number  of  operations  which  must  be  performed  at 
each  time  step  is  proportional  to  the  number  of  '-elLs. 
each  Iteration  requires  less  computer  tune.  Thus,  over¬ 
all  turnaround  time  is  greatly  reduced. 

Futm'e  Plans 

The  method  will  be  applied  to  time-dependent  prublems 
including  those  containing  significant  geometric  com¬ 
plexity.  Evaluation  of  the  performance  of  the  technique 
for  such  problems  is  essential  before  the  technique  will 
be  considered  mature. 

Publications 

Aftosmis.  .M-.T  .  "An  I  pwind  .Method  for  Simulation  of 
Viscous  Flow  on  Adaptively  Refined  Meshes. "  Accepted 
for  publication  by  the  AIAA  Journal.  Mar.  Tjyo. 
Aftosmis.  M..J..  ’’Viscous  Flow  Simulation  Using 
and  Upwind  Method  for  Hexahedral  Based  Adaptive 
Meshes.”  AIAA  Paper  03-0772  Jan,  1993 
Aftosmis.  ”A  Second-Order  T\  D  Method  for 

the  Solution  of  the  3D  Euler  and  Navier-Stokes  Equa¬ 
tions  on  Adaptively  Refined  Meshes."  Proceeding.^  of  tin 
Jjih  InUmaiiOiial  Conference  on  yumerical  Method.'  in 
fluid  Mechanics.  Rome.  Italy,  .]uly  1992- 
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COBALT  -  An  Unstructured  Flow  Solver 

Wiliam  Z  Strang 

Interdisciplinary  and  Applied  Section 


Research  Objective 

To  simplify  computational  fluid  dynamics  (CFD)  anal¬ 
yses  of  air  vehicles  of  arbitrary  geometrical  complexity: 
to  lessen  the  „ime  required  to  generate  such  predictions 
by  a  factor  of  o;  to  reduce  the  CFD  skill  level  required 
of  the  ngineers  conducting  the  analyses. 

A  r  ^roach 

Unstructured  grids  are  one  of  the  very  few  methods,  if 
not  the  only  method  offering  the  potential  to  achieve 
the  above  research  objectives  In  theory,  a  user  gener¬ 
ates  only  the  bounding  surfaces  of  the  problem  and  the 
unstructured  mesh  generator  automatically  tessellates 
the  domain  interior.  Unstructured  grids,  however,  re¬ 
quire  flow  solvers' that  are  capable  of  handling  general 
connectivity  issues  Thus,  an  in-house  project  to  de¬ 
velop  an  unstructured  flow  solver.  COBALT,  began  in 
February  IPPO. 

Accomplishments 

The  fundamental  algorithm  of  COB.ALT  is  the  spatially 
and  temporally  first-order  accurate,  non-linear  schemr 
due  to  Godunov.  The  extension  to  second-order  accu¬ 
racy  in  space  is  due  to  van  Leer  and  Collela. 


Second-  and  third-  order  accuracy  in  time  is  achieved 
with  classical  Runge-Kutta  integration  schemes  Mono¬ 
tonicity,  required  for  stability  of  the  spatially  second- 
order  scheme,  is  based  on  the  works  of  Collela.  Goor- 
jian  and  Obayashi,  Slip  wall,  adiabatic  and  isothermal 
nchsiip  wail.  farfieJd.  inflow  and  ouffiovv  boundary  condi¬ 
tions  are  .available.  Considerable  effort  has  been  devoted 
to  these  boundary  conditions  to  make  the  code  as  robust 
and  as  accurate  as  possible  Construction  of  the  viscous 
terms  follows  MacCormack's  vvork.  The  Baldwin  and 
Barth  turbulence  model  is  used  to  simulate  the  effects  of 
fine  scale  turbulence  Two-dimensional,  axi-syrnmetnc 
and  three-dimensional  problems  may  be  treated.  The 
sole  requirement  of  the  mesh  is  that  it  be  compo.sed  of 
convex  cells.  Lastly,  with  the  exception  of  a  few  pre- 
and  post-processing  tasks.  COBALT  is  fully  vectorized 
and  tvpicallv  operates  in  excess  of  100  MFLOPS  a  Crav 
X-MP/-216  ’ 

Future  Plans 

Additional  validation  of  the  viscous  and  turbulence 
model  routines,  documentation,  third-order  spatial  ac¬ 
curacy.  parallel  implementation,  additional  turbulence 
models  and  grid  adaptation  are  planned 
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Unstructured  Approach  to  the  Design  of  Multiple- Element  Airfoils 

Don  W.  Kinsey 

Interdisciplinary  and  Applied  CFD  Section 


Research  Objective 

Development  of  efficient  high-lift  airfoils  may  lead  to 
reduced  weight  and  complexity  over  current  high  lift 
systems  with  attendant  reduced  cost  and  improved  re¬ 
liability.  The  research  described  here  combines  the  rel¬ 
atively  new  unstructured  grid  computational  fluid  dy¬ 
namics  (CFD)  technique  with  an  existing  inverse  design 
procedure.  The  results  will,  hopefully,  be  a  new  tech¬ 
nique  to  design  and  analyze  high-lift  airfoils  with  the  im¬ 
proved  accuracy  of  an  Euler  solver,  but  significantly  less 
work  for  the  designer  than  most  advanced  CFD  meth¬ 
ods. 

Approach 

The  CFD  solver  vvas  an  implicit  Van  Leer  split  flux  Eu¬ 
ler  scheme  with  point  Gauss-Seidei  relaxation  for  two- 
dimensional  triangular  grids.  A  two-dimensional  un¬ 
structured  grid  generator  based  on  the  Delaunay  tn- 
angulation  process  was  used  to  create  the  grid.  Grid 
movement,  needed  to  allow  the  geometry  to  change  with 
the  inverse  design  procedur-^.  was  provided  by  a  "spring 
analog}'”  process.  The  inverse  design  program  is  a  vari¬ 
ation  of  the  Modified  Garabedian-McFadden  (MGM) 
residiiai-correctiou  procedure.  This  procedure  requires 
a  target  velocity  distribution  and  an  initial  airfoil  ge¬ 
ometry.  Then,  through  a  senes  of  iterations,  the  airfoil 
geometry  is  adjusted  based  on  the  difference  between 
the  current  and  target  veiocit}'  distributions  The  con¬ 
verged  solution  provides  a  geometry  that  produces  the 
target  veiocit.}'  distribution. 


Accomplishments 

The  program  was  developed  and  validated  on  several 
single-element  airfoils.  The  first  multiple-element  airfoil 
test  case  was  a  Whitcomb  slotted  supercritical  airfoil 
with  a  trailing-edge  flap  set  at  20'.  The  target  veloc¬ 
ity  distribution  w-as  obtained  by  direct  Euler  solution  of 
the  known  Whitcomb  geometry.  The  initial  airfoil  ge¬ 
ometry  was  two  symmetric  (N.-VCA  0014)  airfoils  with 
chord  length  and  trailing-edge  locations  of  the  target  ge¬ 
ometries.  The  picture  below  shovvs  the  initial  and  final 
geometries  and  pressure  distributions.  A  three-element 
airfoil  test  case  which  had  a  30“  leading-edge  flap  and  .a 
30“  trailing-edge  flap  was  attempted,  but  proved  to  be 
too  difficult  a  case  for  the  current  method 

Fixture  Plans 

-Ai  three-dimentional  version  of  the  inverse  design  proce¬ 
dure  will  be  developed  The  same  MGM  method  and 
spring  analogy  will  be  used,  however,  a  structured  finite- 
volume  Euler  solver  based  on  Jameson’s  algorithm  wfll 
be  used  for  the  flow  solver.  This  program  will  r-'i.p-iire 
a  two-block  structured  grid  that  is  aligned  streamwi-'e 
over  the  %ving.  A  simple  wing-body  test  case  will  be 
used  to  develop  and  validate  the  code. 

Publications 

Kinsey.  D  W  .  and  Jolly.  B.A..”An  I.’nstructured  .\p- 
proach  to  the  Design  of  Multiple-Element  Airfoils." 
-•MAA  Paper  92-2709.  June  1992. 
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Three-Dimensional  Unstructured  Post-Processing 

Lt  Phillip  M.  Fite 

Interdisciplinary  and  Applied  CFD  Section 


Research  Objective 

Following  grid-generation  and  flow-soluiion.  the  third 
critical  stage  to  the  CFD  solution  process  involves 
post-processing,  or  the  visualization  of  calculated  data 
Post-processing  provides  scientists  and  engineers  with  a 
graphics-based  analysis  tool  for  presenting  CFD  data  as 
well  as  verifying  the  grid-generation  and  flow-solver  ap¬ 
plications.  Due  to  their  simpler  generation  and  quicker 
adaptation  characteristics,  unstructured  grids  are  be¬ 
coming  the  choice  in  the  CFD  community  as  the  need 
to  analyze  more  complex  geometries  increases  To¬ 
day  there  are  many  post-processors  for  structured  grids 
available,  yet  few  meet  the  current  demands  for  visual¬ 
izing  unstructured  grid  data. 

Approach 

An  unstructured  post-processing  capability  is  being  de¬ 
veloped  under  contract  to  Mgyan.  Incorporated  (Dr 
Paresh  Parikh.  Principal  Investigator),  Originally  devel¬ 
oped  for  NASA,  \'igyan's  code  VPLOT3D  was  the  only 
available  post-processor  capable  of  visualizing  unstruc¬ 
tured  :}-D  grids.  Original  versions  of  this  code  were  lim¬ 
ited  to  processing  unstructured  grid  sizes  up  to  IdO.OOO 
grid  points.  To  allow  the  input  of  larger  grids.  Vigyan 
is  modifying  VPLOT3D  by  increasing  the  code  s  mem¬ 
ory  management  capability  and  incorporating  memory 
sharing  techniques  To  improve  memory  management, 
structured  grids  are  handled  separately  from  unstruc¬ 
tured  grids.  Separating  grid  types  reduces  unnecessary 
data  associated  with  structured  meshes,  a.. a  thus  allows 


larger  unstructured  grid  size  This  lncrease^  available 
ntemory  by  reducing  the  need  to  calculate  connectiv¬ 
ity  data  except  when  needed  Incorporating  memory 
sharing  techniques  such  as  computing  derived  quantities 
(Mach  number  .entropy,  etc  i  'on  the  fl>  '  eliminates  the 
need  to  carry  those  values  in  the  processing  arrays 

Accomplishments 

With  the  current  options  available  VPLOT3D  s  un¬ 
structured  capability  has  improved  memory  allocation 
to  process  grid  sizes  of  up  to  2.50.000  grid  points  In 
addition  to  grid  size  input,  working  versions  of  the  <'ode 
possess  some  unique  features  for  unstructured  process¬ 
ing.  These  include  arbirrary  cutting  planes  isovolurne 
plotting,  and  input  of  multiple  data  formats. 

Future  Plans 

To  meet  future  Air  Force  needs  \'PLOT3D  will  be  mod¬ 
ified  to  incorporate  some  additional  features  These  in¬ 
clude  portability  to  multiple  platforms,  time-dependent 
visual  simulation  i animation'  of  unsteady  flows  and 
moving  objects,  plotting  time  dependent  data  on  a  mov  ¬ 
ing  plane  or  grid,  plotting  experimeiita!  data,  and  store 
separation  analysis 

Publications 

Parikh.  P..  et  al.  A  Package  For  3-D  Fnstructured  Grid 
Generation.  Finite  Element  Flow  Solution  and  Flow 
Field  Visualization."  NASA  Contractor  Report  182090. 
September  1990. 
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TOPDUUG  -  The  Only  Package  for  Design  Using  Unstructured  Grids 

Marvin  C.  Gridley 

Airframe  Propulsion  Weapons  Integration  Section 


Research  Objective 

The  use  of  CFD  in  research  and  design  can  often  b  *  a 
time-consuming  and  somewhat  esoteric  "art”  to  scien¬ 
tists  and  engineers  who  may  not  be  completely  familiar 
with  the  inner  workings  of  grid  generators,  flow  solvers, 
and  post-processors.  Highly  interactive  computer  pro¬ 
grams  which  provide  graphical  feedback  have  proven  to 
be  great  tools  for  use  in  research.  This  effort  focused 
on  the  generation  of  a  computer  program  designed  to 
integrate  se^e^al  of  the  tasks  in  the  CFD  process  into 
an  interactive,  graphical  package  which  was  easily  used 
by  any  engineer  or  scientist,  not  just  by  CFD  code  de¬ 
velopers. 

Approach 

Several  different  computer  codes  were  being  used  to  per¬ 
form  the  tasks  of  grid  generation  and  pre-  ano  post¬ 
processing  within  the  Aeromechanics  Divisi^jn.  These 
codes  ivere  combined  in  TOPDUUG  to  create  a  general 
"cradle-to-grave"  tool  for  use  in  2D  CFD  analysis.  Inter¬ 
active  menus  and  windows  were  added  to  increase  func¬ 
tionality  and  utiiitj'.  An  interactive  on-line  help  system 
w^as  incorporated  Capabilities  of  TOPDUUG  include: 

1.  Fully  interactive  boundary  condition  specification 

2.  Two-dimensional  unstructured  grid  generation, 
including  viscous  grids 


3.  Flow  solver  job  file  generation 

4.  Flow  solution  visualizaU  n 

5.  Time  dependent  flow  s  jtion  animation 

Accomplishments 

The  first  production  version  of  TOPDUUG,  The  Only 
Package  for  Design  Using  Unstructured  Grids,  has  been 
produced.  TOPDUUG  is  a  highly  interactive  2D  un¬ 
structured  grid  generator/pre/post-processor  designed 
for  ease  of  use  and  utility.  TOPDUUG  is  currently  be¬ 
ing  used  in  the  Aeromechanics  Division  for  several  en¬ 
gineering  applications  as  well  as  for  code  development 

Future  Plans 

Future  plan.=  include  modification  of  the  2D  unstruc¬ 
tured  grid  generation  module  to  be  more  robust  and 
user  friendly.  Several  scalar  and  vector  functions  will  be 
added  to  the  post-processing  module  to  increase  func¬ 
tionality.  A  user's  manual  describing  operation  of  the 
softrvare  will  be  written. 

Publications 

Gridley.  M,  C  .  ‘  TOPDUUG  -  The  Only  Package  for  De¬ 
sign  Using  Unstructured  Grids  User's  Manual."  to  be 
published. 


Several  Screens  Showing  TOPDUUG  Functions 
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Numerical  Simulation  of  Turbulent  Cylinder  Juncture  Flowfields 

Donald  P.  Rizzetta 
CFD  Research  Section 


Research  Objective 

This  work  investigates  the  ability  of  an  advanced  two- 
equation  turbulence  model  to  properly  simulate  flow- 
fields  which  exhibit  complex  three-dimensional  separa¬ 
tions.  Steady  high  Reynolds  number  subsonic  and  su¬ 
personic  flows  past  a  cylinder  which  is  mounted  upright 
on  a  flat  plate  are  considered.  These  flowfields  are  char¬ 
acterized  by  formation  of  a  horseshoe  vortex  system  that 
wraps  around  the  cylinder,  and  are  of  fundamental  in¬ 
terest  because  of  the  intricate  fluid  phenomena  which 
occur.  In  addition,  the  configuration  is  representative  of 
a  variety  of  practical  situations  including  wing/fuselage 
and  wing/pylon  intersections,  struts  of  aerospace  vehi¬ 
cles.  submarine  conning  towers,  turbomachinery,  wind 
tunnel  model  supports,  architectural  aerodynamics,  and 
meteorological  and  geological  applications. 

Approach 

Cylinder/plate  juncture  flowfields  are  simulated  nu¬ 
merically  by  solution  of  the  time-dependent  three- 
dimensional  compre.ssible  mass-averaged  Navier-Stokes 
equations,  which  are  integrated  in  time  to  achieve  the 
steady  state.  Effects  of  fine  scale  turbulence  are  repre¬ 
sented  by  a  two-equationlF  —  ()  closure  model  which  in¬ 
cludes  a  generalized  formulation  and  low- Reynolds  num¬ 


ber  terms.  For  supersonic  flow,  the  turbulence  equation.s 
incorporate  a  compressibility  correction. 

Accomplishments 

An  existing  vectorized  three-dimensional  implicit  com¬ 
puter  code  was  modified  in  order  to  accommodate  so¬ 
lution  of  the  turbulence  equations,  and  the  code  was 
validated.  Both  subsonic  and  supersonic  cylinder/plate 
juncture  flowfields  were  generated  numerically.  All  com¬ 
monly  observed  physical  features  of  the  flows  were  repro¬ 
duced  by  the  calculations.  A  grid  mesh  step-size  study 
was  performed  in  order  to  assess  resolution  requirements 
of  the  solutions,  and  the  effect  of  the  compressibility- 
correction  was  examined  Results  of  the  computations 
have  been  compared  to  experimental  data  in  terms  of 
static  and  total  pressure  coefficients,  velocity  profiles, 
and  surface  limiting  streamline  patterns. 

Future  Plans 

The  turbulence  model  equations  will  be  used  to  simulate 
the  flow  over  a  delta  wing  performing  a  roll  maneuver, 
as  well  as  the  explosive  rupture  of  an  aircraft  forebody. 

Publications 

Rizzetta.  D  P  “Numerical  Simulation  of  Turbulent 
Cvlinder  Juncture  Flowfields."  AlAA  Paper  93-3038, 
July  1993. 


Turbulence  Modeling  for  Thrust  Reverser  Flow  Prediction  Methods 


Kenneth  E.  VVurtzler 

Interdisciplinary  and  Applied  CED  Section 
Research  Objective 

Thrust  vectoring  is  one  means  to  improving  the  ma- 
neu%erabiiity  and  take-off  and  landing  performance  of 
aircraft.  However,  thrust  vectoring,  especially  in  close 
proximity  to  the  ground,  can  significantly  increase  the 
complexity  of  the  flowfield  around  the  aircraft.  Prob¬ 
lems  that  can  result  include  hot  gas  ingestion,  reduced 
lift,  and  foreign  object  damage  to  engines.  To  achieve 
good  aircraft  designs,  it  is  necessary  to  have  an  accurate 
now  prediction  method  capable  of  resolving  turbulence. 
This  work  addresses  the  issue  of  turbulence  modeling 
for  the  inclined  impinging  jets  associated  with  the  use 
of  thurst  reversing  in  ground  effect. 

Approach 

This  effort  combined  numerical  research  with  experi¬ 
mental  results  in  order  to  produce  a  validated  turbu¬ 
lence  model.  The  cases  studied  were  an  impinging  jet 
issuing  from  a  nozzle  normal  to  the  ground  plane  and 
one  issuing  from  a  nozzle  inclined  45'’  into  the  cross- 
flow.  Nielsen  Engineering  and  Research.  Inc.  (Dr. 
Robert  Childs.  Principal  Investigator)  performed  large 
eddy  simulations  (LES)  to  obtain  mean  flow  and  tur¬ 
bulence  data.  Basic  physical  experiments,  performed  at 
Stanford  University,  provided  additional  flow  data. 


This  data  was  analyzed  to  identify  the  dominant  phys¬ 
ical  processes  which  govern  turbulence  in  these  flows. 
Modifications  to  the  k-e  turbulence  model  were  based 
on  these  findings  The  new  coefficients  in  the  resulting 
modified  turbulence  model  were  determined  via  numer¬ 
ical  optimization  for  benchmark  jet  flows. 

Accomplishments 

Streamline  curvature  and  vortex  stretching  were  identi¬ 
fied  as  the  two  most  significant  influences  on  turbulence. 
Modifications  which  account  for  these  two  features  were 
added  to  the  standard  k-e  turbulence  model.  In  the 
cases  considered,  the  modified  model  consistently  im¬ 
proved  the  accuracy  for  the  computed  location  of  the 
ground  vortex.  Errors  in  the  computed  vortex  location 
were  as  large  as  50%  in  calculations  with  the  standard 
turbulence  model  and  could  be  of  either  sign.  The  mod¬ 
ified  model  reduced  these  errors  to  the  range  of  13-18%. 

Publications 

WL-TR-92-3I21.  "Turbulence  Modeling  for  Thrust  Re¬ 
verser  Flow  Field  Prediction  Methods”.  December  1992. 

Childs.  R.  E.  and  Rodman.  L  C  .  "Turbulence  Modeling 
for  Impinging  Jet  Flows.”  AIAA-92-2672. 
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The  Effect  of  Leading-Edge  Cross-Sectional  Geometry  on  Vortex  Flow 
Aerodynamics 


T.  Sean  Tavares  -  NRC  Associate 
CFD  Research  Section 

Research  Objective 

The  objective  of  this  research  is  the  detailed  calculation 
of  viscous  flow  features  which  govern  vortex  separation 
from  the  leading  edges  of  highly  swept  wings  at  high 
Reynolds  numbers.  The  extent  to  which  leading-edge 
vortex  separation  occurs  can  have  a  strong  effect  on  lift, 
drag,  and  pitching  moment,  and  depends  in  part  on  the 
leading-edge  cross-sectional  geometry.  For  wings  with 
sharp  leading  edges,  the  point  of  primary  vortex  sepa¬ 
ration  is  fixed  at  the  leading  edge.  However,  for  wings 
having  rounded  leading  edges,  this  is  not  necessarily  the 
case,  and  the  flow  may  remain  attached  along  part  or 
all  of  the  leading  edge  depending  on  angle  of  attack. 
Reynolds  number,  and  Mach  number. 

Approach 

The  approach  taken  is  the  computational  solution  of  the 
3-D  Reynolds- .Averaged  .Navier-Stokes  equations.  They 
are  solved  using  the  implicit,  approximately  factored 
Beam-Warming  algorithm  as  implemented  within  an  in- 
house  computer  code  (FDL3DI).  Turbulence  closure  is 
achieved  vvith  the  Baldwiu-Lomax  algebraic  turbulence 
model,  which  has  been  modified  to  account  for  the  ef¬ 
fects  of  large  scale  cro.ssflow  separation.  A  special  chal¬ 
lenge  of  this  problem  is  to  generate  a  grid  with  sufficient 
resolution  to  resolve  the  thin  boundary  layers  character¬ 
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istic  of  high  Reynolds  number  flows  and  to  maintain 
good  grid  quality  in  the  rest  of  the  domain,  without 
using  an  e.xcessive  number  of  grid  points. 
Accomplishments 

The  figure  illustrates  a  calculation  of  the  flow  about  a 
clipped  delta  wing  with  65°  leading-edge  sweep  and  a 
sharp  leading  edge.  The  flow  conditions  are  10°  an¬ 
gle  of  attack  at  Mach  0.4  and  a  Reynolds  number  of 
9  X  10°.  Good  agreement  between  the  calculated  and 
measured  spaiiwise  pressure  distributions,  as  shown  for 
the  station  at  60%  of  the  root  chord  is  indicative  of 
the  successful  calculation  of  both  the  primary  and  sec¬ 
ondary  vortex  structures.  These  calculations  represent 
the  highest  Reynolds  number  turbulent  flow  cases  run 
to  date  with  the  in-house  Navier-Stokes  code,  hi  addi¬ 
tion.  comparisons  of  preliminary  solutic-.r  obtained  on 
coarse  grids  for  both  sharp  and  round-edged  wings,  r**- 
produce  observed  dependences  of  vortex  flow  patterns 
on  leading-edge  geometry. 

Future  Plans 

Calculations  are  underway  on  the  65°  clipped  delta  with 
rounded  leading  edges,  for  which  experimental  measure¬ 
ments  show  separation  along  only  part  of  the  leading 
edge  at  10°  angle  of  attack.  Calculations  are  planned 
for  both  geometries  on  refined  computational  meshes 
having  1.6  million  grid  points. 


Prvnuw  OMibuion  at  sm  Chant 
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Vortex  Breakdown  Above  a  Pitching  Delta  Wing 

Mi^el  R.  Visbal 
CFD  Research  Section 


Research  Objective 

V'ortex  breakdown  represents  one  of  the  unresolved  is¬ 
sues  in  high-angle-of-attack  aerodynamics.  Further  un¬ 
derstanding  of  this  complex  fluid  dynamics  phenomenon 
is  required  for  enhancing  aircraft  maneuverability  and 
for  alleviating  undesirable  fluid/structure  interactions 
such  as  tail  buffet.  The  purpose  of  this  research  is  to 
study  the  transient  behavior  and  3-D  flow  structure  of 
the  leading-edge  vortex  breakdown  above  a  delta  wing 
which  is  pitched  to  a  high  angle  of  attack. 

Approacli 

The  flows  are  simulated  by  solving  the  unsteady  three- 
dimensional  full  Navier- Stokes  equations  expressed  in 
strong  conservation  law  form.  A  general  time-dependent 
coordinate  transformation  is  incorporated  to  treat  the 
prescribed  wing  motion.  The  discretized  equations  are 
solved  employing  the  implicit  approximate-factorization 
Beam- Warming  algorithm.  The  accuracy  and  stability 
of  the  scheme  are  enhanced  by  means  of  a  subiteration 
procedure. 

Accomplishments 

Calculations  were  performed  for  a  75°  sweep  delta  wing 
pitching  at  a  constant  non-dimensional  pitch  rate  Q"*"  = 
0.3.  from  an  initial  angle  of  attack  a,  =  25°  to  a  final 
angle  a;  —  50° .  The  pitch  a.xis  was  located  at  the  wing 
trailing  edge  .An  assessment  of  the  effects  of  numerical 
resolution  and  the  favorable  comparison  with  available 
experimental  data  indicate  the  computational  approach 
captures  the  basic  dynamics  of  the  transient  breakdown. 


The  angular  delay  and  onset  of  breakdown  are  found 
to  be  strongly  linked  to  the  pressure  gradient  prevailing 
along  the  vortex  axis  A  description  of  the  3-D  instan¬ 
taneous  structure  of  vortex  breakdown  is  provided  for 
the  first  time  using  critical-point  theory.  The  region  of 
reversed  flow  in  the  vortex  core  is  associated  with  the 
appearance  of  pairs  of  opposite  3-D  spiral/saddle  crit¬ 
ical  points.  During  its  early  stages,  the  vortex  break¬ 
down  is  fairly  axisymmetric.  However,  as  breakdown 
proceeds  upstream,  asymmetric  effects  become  increas¬ 
ingly  important  and  lead  eventually  to  the  formation  of 
a  breakdown  bubble.  This  bubble  structure  is  open  and 
contains  within  itself  a  pair  of  stagnation  points  which 
are  diametrically  opposed  and  which  rotate  in  the  same 
sense  as  the  base  flow.  .4  representation  of  the  break¬ 
down  bubble  using  an  iso-surface  of  constant  total  pres¬ 
sure  shows  a  great  deal  of  resemblance  to  experimental 
flow  visualizations  of  axisymmetric  vortex  breakdown  in 
a  tube. 

Future  Plans 

Further  investigation  of  the  relation  between  the  dif¬ 
ferent  representations  of  the  breakdown  region  (i.e.  in¬ 
stantaneous  flow  structure,  streakline  visualizations  and 
mean-flow  measurements)  will  be  pursued  for  unsteady 
breakdown  above  a  stationary  delta  wing.  Investigation 
of  the  parametric  effects  of  pitch  rate  and  pitch-axis  lo¬ 
cation.  as  well  as.  compressibility  effects  is  in  progress. 
Publications 

Visbal.  M.  R..  “Structure  of  Vortex  Breakdown  on  a 
Pitching  Delta  Wing.”  AIAA  Paper  93-0434.  January 
1993. 


Computed  bubble-type  breakdown  above  pitching  delta  wing 
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Numerical  Simulation  of  Delta-Wing  Roll 

Raymond  E.  Gordnier 
CFD  Research  Section 


Research  Objective 

The  combat  requirements  for  modern  fighter  aircraft  ne¬ 
cessitate  that  they  be  highly  maneuverable  and  able  to 
operate  at  high  angles-of-attack.  These  requirements 
lead  to  an  interest  in  understanding  the  aerodynam¬ 
ics  of  delta-wing  maneuvers,  including  high-rate  rolling 
motion.  The  focus  of  the  present  research  effort  is  the 
numerical  simulation  of  a  constant  roll-rate,  “roll-and- 
hold”  maneuver.  The  main  objective  is  to  describe  the 
dynamical  behavior  of  the  vortex  position  and  strength, 
as  well  as  the  corresponding  effect  on  surface  pressure, 
lift  and  roll  moment. 

Approach 

The  governing  equations  solved  are  the  unsteady,  three- 
dimensional,  full  Navier-Stokes  equations  in  strong  con¬ 
servation  form.  A  general  time-dependent  coordinate 
transformation  is  introduced  to  allow  for  the  body  mo¬ 
tion.  The  equations  are  solved  using  the  implicit,  ap¬ 
proximately  factored  Beam- Warming  algorithm.  Both  a 
full  block  tridiagonal  and  diagonalized  inversion  schemes 
are  available.  A  subiteration  procedure  is  also  incorpo¬ 
rated  to  improve  the  accuracy  and  robustness  of  the 
algorithm. 

Accomplishments 

Two  fixed-roll  angle  cases  (0  =  0°  and  0  =  45®) 
were  computed  to  demonstrate  the  capabilities  of  the 


numerical  scheme  before  computing  the  dynamic  roll 
maneuver.  The  results  showed  excellent  qualitative 
and  quantitative  agreement  with  experimental  measure¬ 
ments.  Then  a  constant  rc^l  rate,  4*  =  0.1325.  “roll-and- 
hold”  maneuver  from  ©  =  0®  to  0  =  45®  was  computed. 
The  dynamics  of  the  vortices  during  the  roll  maneuver 
have  been  described.  A  lag  in  the  upward  leading-edge 
vortex  similar  to  that  observed  for  delta  wing  roA  was 
noted.  The  effect  of  the  dynamic  motion  on  the  vortex 
strengths  has  also  been  described-  A  simple,  quasi-static 
explanation  of  this  vortex  behavior  based  on  effective 
angle-of-attack  and  sideslip  angle  during  roll  has  been 
proposed. 

Future  Plans 

Current  work  on  this  project  consists  of  the  investigation 
of  parametric  effects  on  the  aerodynamics  of  the  roll  ma¬ 
neuver.  The  effects  of  roll  rate.  Reynolds  number,  and 
asymmetric  initial  conditions  (i.e..  starting  from  nonzero 
roll  angles)  are  being  studied.  Calculations  of  turbulent 
flows  over  delta  wings  using  a  two  equation  k-c  turbu¬ 
lence  model  are  also  planned.  The  ultimate  objective 
of  this  work  is  to  be  able  to  numerically  simulate  delta- 
wing  rock  both  with  and  without  vortex  breakdown. 
Publications 

Gordnier.  R.  E.  and  Visbal.  M  R..  “Numerical  simula¬ 
tion  of  Delta- Wing  Roll."  AIAA-93-0554.  January.  1993. 
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Nonequilibrium  Hypersonic  Flowfield  Prediction  for 
Flow  Past  Blunt  Bodies 


M^war  .I>)syiila 
(  i  Hi'si'arcli  Si'ctio!) 

R.<!s<'ai'ch  Ohjt!<!tiv(‘ 

IV)  p.-rl'orm  iiuiiKTu  a!  niod.  llmp;  of  tlx-  ilowfixl.l  using 
Iho  \avi-'r-Stok('s  i'<(uat  i.)tis  with  ililili'  ralx  rlxanical  ki- 
iixtics  applioahh'  to  lugli  t('tiip<‘iatur('  air  (or  liypi'rsDitir 
(lows  ill  eh  t'niical  and  tln'nnal  noinapnliltriiini. 
Approach 

111!'  axisyinmct  lie  Navii'r-Stok‘'s  Ivpiations.  •'>  spiv-ics 
coniimiily  ('<|iialioii.s,  ;i  viliratioiial  cni-rgy  cipialioiis, 
and  a  total  energy  e(|ualion  are  s<jlved  in  a  fully  cou- 
|ded  explicit  met  hod  using  t  he  Hoe  I'liix-niirerenc,'  Split 
Scheme, 

Accouiplisiiiiimits 

I  he  extension  of  till-  Hoe  scheme  to  incinde  the  finile 
rate  chenucal  kineiic  eipiations  follows  the  approach 
of  l  irossman  and  Cinnella.  Second  ordc-c.' accuracy  is 
achieverl  with  the  Mf’SCh  a|)proach  Tlu-  full  Navier- 
Stoke.s  equations  are  solved  with  (inile  rate  chemistry 
for  the  (low  past  an  axi-synunetric  blunt  body  at  zero 
incidence  at  several  Mach  numbers  ranging  from  Mach 
10  to  18,  at  altitudes  of  22  km  and  37  km.  The  wall 
temperature  is  set  at  555  K  and  the  maximum  shock 
temperature  is  around  15.000  K.  Both  non-catalytic  and 
fully  catalytic  wall  boundary  conditions  for  the  species 
are  utilized. 


rile  air  simulated  m  tin-  modetimg  is  composed  of 
aioimc  oxygen  and  nitrogen,  moleciitar  oxygen  and  ni¬ 
trogen  and  nitric  oxide  flie  cf.emical  source  terms 
are  derived  I'roiii  the  Law  of  .Mass  .■Action  Fhe  im¬ 
portance  of  a  mull i-ieiiiperat lire  model  for  the  proper 
Irealnienl  of  the  viliral  lotial  relaxation  in  the  predn 
lion  of  the  surface  heat  transfer  is  shown.  I’lie  variation 
in  compiitetl  shock -standoff  distance  snbslantiate.s  tlial 
the  .Mach  number  indepi-ndence  principle  a|)plicah|e  to 
ideal  gasi's  does  not  apply  to  dis.sociating  flows.  I  he 
computed  stagnation  point  heat  flux  and  tlie  surface 
heat  transfer  show  excellent  comparison  with  the  ex¬ 
perimental  ilaia  obtained  from  recent  exiierimeiits  con¬ 
ducted  in  a  shock  tunnel  and  from  a  set  of  earlier  shock 
tube  experiments. 

Fiitmxi  Plans 

l-'dfects  of  ionization  and  radiation  become  important  at 
higher  temperatures  and  they  will  he  studied  for  proper 
prediction  of  the  flowlield  and  surface  heat  transfer  on 
the  How  past  blunt  l>odies 

Publications 

Josyula,  Eswar  and  Shang,  Joseph  “(,’ompula- 
tion  of  Nonequilibrium  Hypersonic  Flowfields  Past 
Hemisphere-Cylinders.”  Journal  of  Thermophysics  and 
Heat  'Transfer,  to  be  pnbiisiiea  in  July  1993. 
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Stability  Analysis  of  a  Combined  Couette-Poiseuille,  Two-Fluid  Flow 


Lt  John  J.  NVlson 

Interdisciplinary  and  Applied  CFD  Section 
Research  Objective 

In  this  work,  the  complete  theory  is  developed  for  the 
linear  stability  analysis,  energy  analysis  and  weakly  non¬ 
linear  analysis  for  a  combined  Couette-Poiseuille,  two 
fluid  flow.  The  linear  stability  analysis  determines  the 
conditions  at  which  the  basic,  laminar  flow  will  become 
unstable  and  transition  occurs.  The  energy  analysis  de¬ 
termines  the  mechanisms  which  cause  the  instability. 
In  multi-fluid  incompressible  flow,  instability  is  influ¬ 
enced  by  five  mechanisms:  viscous  dissipation,  Reynolds 
stress,  interfacial  friction,  density  stratification  and  sur¬ 
face  tension.  Weakly  nonlinear  analysis  determines  if 
standing  wave  solutions  are  possible  in  the  flow  ensu¬ 
ing  after  Iransistion,  and  if  so,  determines  the  transient 
amplitude  of  the  waves. 

Approach 

The  linear  stability  analysis  is  accomplished  by  linearly 
perturbing  the  continuity,  Navier-Stokes  and  interfacial 
equations.  A  separable  solution  is  assumed,  which  leads 
to  a  complex  eigensystem  problem,  which  determines 
the  growth  rate  of  infinitesimal  disturbances  in  the  flow. 
Full  viscous  amd  surface  tension  effects  are  included. 
The  results  of  the  linear  analysis  are  used  to  derive  a 
enertry  ane'ysis.  The  weakly  nonlinear  analysis  is  ac¬ 
complished  by  expanding  in  a  small  parameter  to  third 


k 

(a) 


order  about  the  basic  flow  Systems  governing  the  di-s- 
toriion  of  the  mean  flow,  the  rise  of  the  second  itarinonic 
and  growth  rate  of  finite  amplitude  waves  are  found  and 
solved. 

Accomplishments 

This  work  completely  derives  a  theory  and  provides  test 
case  solutions  for  transition  prediction  coding  of  incom¬ 
pressible  flows.  The  flow  situations  studied  in  the  work 
are  the  standard  test  cases  for  linear  stability  code  accu¬ 
racy.  Our  linear  stability  analysis  is  in  agreement  with 
the  standard  answers.  Going  beyond  linear  analysis, 
this  work  is  the  first  to  identify  the  mechanisms  leading 
to  the  instability.  Specific  questions  relating  to  normal¬ 
izations  are  also  addressed  for  the  first  time. 

Future  Plans 

The  theory  developed  is  being  employed  in  an  incom¬ 
pressible  Navier-Stokes  code  and  used  to  study  a  prob¬ 
lem  of  current  interest.  Also  planned  is  the  development 
of  a  stability  theory  for  compressible  flows.  This  theory 
will  be  used  to  give  transition  prediction  capability  to 
our  in-house  compressible  Navier-Stokes  code 


(a)  Linear  stability  analysis  of  two  fluids  in  a  Couette  flow.  R  represents  the  Reynolds  number  and  k  the  wave 
number,  (b)  F.nergy  analysis  for  the  flow  of  (a)  at  R  =  600  The  term  "body  term  "  represents  viscous  dissipation 
plus  Reynolds  stress. 
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Computational  Aerodynamic  Analysis  of  a  Decoy  Configuration 

Frank  C.  Witzeman  Jr 
Interdisciplinary  and  Aj^plied  CFD  Section 


Research  Objective 

Aerodynamic  analyses  of  air  vehicles  operating  over  a 
range  of  flight  conditions  can  be  performed  by  a  va- 
riet\  of  methods  including  analytical,  semi-empirical. 
Computational  Fluid  Dynamics  (CFD).  and  experimen¬ 
tal  testing  Limitations  in  both  the  semi-empirical  and 
testing  techniques  for  a  decoy  configuration  lead  to  the 
use  of  CFD  to  fill  gaps  in  the  required  data.  Specif¬ 
ically.  aerodynamic  data  at  transonic,  high  angle-of- 
attack  conditions  were  needed.  Primary  interest  was 
in  obtaining  in-plane  (pitch)  forces  and  moments  due  to 
inertial  effects  (inviscid  flow).  Effects  due  to  roll  angle 
and  small  changes  in  Mach  number  and  angle-of-attack 
were  of  interest,  but  secondary. 

Approach 

The  in-house  Euler  code  MERCURY  was  used  to  obtain 
CFD  solutions  for  the  decoy.  A  structured  grid  consist¬ 
ing  of  -15  blocks  and  approximately  720.000  grid  points 
was  generated  from  the  GRIDGEN  code.  The  block¬ 
ing  scheme  wa.s  beneficial  in  reducing  computer  mem¬ 
ory  requirements  where  the  maximum  memory  required 
depends  only  on  the  largest  block.  Integrated  and  com¬ 
ponent  forces  and  moments  were  compiled  for  each  so¬ 
lution.  The  components  consisted  of  the  nose,  center- 
body.  aft-body,  and  tail.  The  integrated  and  component 


results  were  compared  with  the  semi-empirical  method 
results  and  experimental  data  (integrated  only).  The 
final  solution  matrix  consisted  of  one  roll  angle,  three 
Mach  numbers,  and  several  angles-of-attack. 

Accomplishments 

Euler  solutions  for  the  decoy  at  Mach  numbers  of  0.5 
and  0  7  have  been  obtained  for  several  angles-of-attack 
and  O'"  roll  angle.  The  low-speed  results  compare  well 
to  the  experimental  ditta  up  to  intermediate  angles-of- 
attack  (20°).  whereas  the  Mach-0.7  results  show  discrep¬ 
ancies  beyond  10°  angle-of-attack.  Compressibility  ef¬ 
fects  have  a  major  impact  on  the  aerodynamics  as  noted 
in  the  CFD  results  These  results  also  indicated  that 
the  semi-empirical  method  did  not  accurately  predict 
cenier-of-pressure  location,  especially  for  the  aft-body 
component 

Future  Plans 

Additional  CFD  solutions  at  Mach- 1.2  and  angles-of- 
attack  up  to  90°  are  planned.  These  results  are  lacking 
in  the  experiments.  All  results  will  be  tabulated  for 
the  total  integrated  and  component  forces  and  moments. 
Compressibility  effects  at  transonic  conditions  will  be 
investigated  in  more  detail  with  the  CFD  solutions. 


.\I.icb  Number  (.'ontours  .and  Vi-locity  Vectors  for  Decoy  ( 'onfignr.otion 
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Euler  Analysis  for  Delta  Wing  Design 

Lt  John  S.  Seo 

Interdisciplinary  and  Applied  CFD  Section 


Research  Objective 

At  the  request  of  ASC/XR,  the  feasibility  of  a  14%  thick 
flying  delta  wing  to  sustain  efficient  flight  at  a  design 
speed  of  Mach  0.90  and  lift  coefficient  of  0.60  was  tested. 
CFD  analysis  was  applied  for  this  configuration  as  the 
calculation  of  wave  drag  using  empirical  methods  were 
unreliable.  This  delta  wing  has  a  leading-edge  sweep 
angle  of  53.1  degrees  and  half- wing  planform  area  of 
523.87  square  feet. 

Approach 

Since  the  parameter  of  interest  was  wave  drag,  an  Euler 
solution  was  sufficient  for  this  analysis.  Three  cross- 
sectional  cuts  of  the  wing  were  prepared  in  an  interac¬ 
tive  graphics  grid,  generator.  I3G/VIRGO,  to  create  a 
full  three-dimensional  geometric  surface  definition.  A 
single  C-H  type  grid  with  dimension  of  61x40x70  was 
created  using  another  interative  multiple- block  grid  gen¬ 
erator.  GRIDGEN.  MERCURY,  an  in-house  3-D  Euler 
flow  solver  based  on  Jameson's  algorithm,  was  used  for 
wave  drag  calculations. 


Accomplishments 

Over  50  solutions  with  Mach  number  ranging  from  0.5 
to  0.95  and  angle-of-attack  ranging  from  0  to  15  de¬ 
grees  were  computed  on  two  of  AFMC's  supercomput¬ 
ers:  Wright-Patterson  AFB's  Cray/XMP  .and  Eglin 
AFB's  Cray/YMP.  It  was  found  that  at  Mach  =  0.90, 
and  even  at  a  lower  Mach  of  0.80.  this  14%  thick  delta 
wing  will  not  achieve  the  desired  C;  of  0.60.  A  strong 
shock  slightly  past  the  leading  edge  and  reversed  flow 
near  the  wing  tip  contribute  to  the  Ci  loss 
Significance 

E.xisting  CFD  tools  allow  VVL/FIMC  to  perform  an  Eu¬ 
ler  study  for  this  relatively  simple  configuration  with 
ease.  This  particular  study  lasted  5  weeks  from  begin¬ 
ning  to  end 
Presentations 

The  results  of  this  analysis  were  presented  at  the 
Dayton-Cincinnati  AI.4A  mini-symposium,  25  Mar  93. 


MERCERY  solution  for  delta  vving  at  Mach  0  80  and  8  degrees  angle  of  attack 


I 
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Numerical  Simulation  of  a  Jet  Expulsion  from  an  Aircraft  Forebody, 
With  Comparison  to  Experiment 

Capt  James  Mu'idy  -  CFD  Research  Section 


Research  Objective 

This  research  investigates  the  effect  a  sonic,  underex- 
panded  jet  issuing  from  a  generic  aircraft  forebody,  has 
on  the  airflow  around  the  aircraft,  and  on  the  resulting 
aerodynamic  forces. 

Approach 

An  implicit  Beam- Warming  finite  difference  scheme  is 
being  used  to  solve  the  Navier-Stokes  equations  around 
the  aircraft  forebody.  A  k  —  (  turbulence  model  is  used 
to  account  for  the  effects  of  turbulence.  The  result¬ 
ing  numerical  solutions  are  compared  with  pressure  data 
and  flow  visualization  video  and  still  photography  taken 
from  wind-tunnel  tests.  Pressure  data  is  in  two  forr  s. 
the  first  being  pressure  tap  data,  the  second  being  data 
taken  from  a  new  pressure  sensitive  paint  technique. 
Accomplishments 

The  three-dimensional  Navier-Stokes  equations,  vvith  a 
k  —  €  turbulence  model  have  been  solved  around  the 
generic  aircraft  forebody  on  a  108x81? x43  grid  (408.6T2 


grid  points),  without  the  jet  blowing.  Agreement  be¬ 
tween  the  numerical  solution  and  experimental  data  is 
excellent.  The  solution  for  the  jet-on  case  is  still  only 
partially  converged,  though  indications  are  that  the  nu¬ 
merical  solution  will  agree  closely  with  the  experimental 
data 

Future  Plans 

The  jet-off  case  wdll  be  refined  using  a  finer  grid  (ap¬ 
proximately  1,25  million  grid  points)  to  provide  valida¬ 
tion  for  the  CFD  scheme.  Once  this  is  done,  the  generic 
aircraft  forebody  will  be  modeled  to  conform  to  typical 
flight  characteristics  of  a  jumbo-jet.  The  free-flight  case 
will  mode!  flight  at  Mach  0.7.5.  at  an  altitude  of  30.000 
feet. 

Publications 

Melville,  Reid  "Numerical  Simulation  of  a  Sonic.  Un¬ 
derexpanded  .Jet  from  an  Aircraft  Forebody"  FA.U  Air¬ 
craft  Hardening  and  Survivability  Symposium.  August 
1992.  Atlantic  City  N'J 


Streamlines  Around  a  Jet  Issuing  from  a  Fuselage 
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Numerical  Analysis  of  a  Chined  Forebody  with  Asymmetric  Strakes 


Kenneth  E.  Wurlzler 

Interdisciplinary  and  Applied  CFD  Section 
Research  Objective 

Vortex  forebody  control  (VFC)  is  increasingly  being 
viewed  as  a  way  to  provide  greater  directional  control 
of  fighter  aircraft  at  high  angle  of  attack.  The  major 
objective  of  this  work  to  examine  the  effectivenes.s  of 
various  strake  sizes  on  the  yawing  and  rolling  moments 
on  the  Fighter  Lift  and  Control  (FLAC)  forebody  alone. 
A  secondary  purpose  is  to  validate  the  use  of  an  Euler 
code  for  the  prediction  of  flow  over  a  chined  forebody  at 
high  angles  of  attack. 

Approach 

The  geometry  utilized  was  the  full  FLAC  forebody  with 
an  attached  sting  extending  behind  the  final  fuselage 
cut.  Different  sized  strakes  were  attached  to  the  left 
side  at  the  apex  of  the  nose.  A  1.7  million  point,  multi- 
block  grid  w<is  constructed  using  I3G\TRGO  In  order 
to  simplify  the  grid  generation  and  to  provide  for  quick 
substitution  of  different  strake  sizes,  the  ability  to  create 
various  strake  sizes  was  built  into  the  grid.  Grid  spacing 
normal  and  tangential  to  the  surface  was  kept  equal  so 
as  to  produce  a  dense,  equally  spaced  gridded  domain 
capable  of  capturing  the  core  of  the  vortex.  The  flow- 
past  the  forebody  was  computed  by  using  an  in-house 
Euler  code.  MERCURY.  The  full  'forebody  was  mod¬ 
eled  so  as  to  capture  any  asymmetry.  Convergence  on 
the  forebody  alone  case  was  reached  in  about  2.000  it¬ 
erations.  Each  strake  case  needed  approximately  3.000 
more  iterations 


Convergence  was  determined  by  a  drop  of  four  orders  of 
magnitude  of  the  residuals  and  leveling  off  of  the  body 
forces. 

Accomplishments 

Solutions  were  obtained  at  15.  25.  35.  and  45°  angle  of 
attack  at  a  Mach  number  of  018  Initial  calculations 
on  the  clean  forebody  resulted  in  no  side  force  being 
produced,  which  is  theoretically  correct  for  Euler  codes. 
The  results  showed  that  a  short,  narrow  strake  produced 
a  yawing  force  that  acts  m  the  direction  of  the  strake 
side.  Wider  strakes  pushed  the  vortex  away  from  the 
forebody,  creating  a  yawing  force  that  acts  in  the  direc¬ 
tion  of  the  other  side.  The  strong  vortex  that  separates 
from  the  chined  forebody  is  captured  vvell  by  the  Euler 
code 

Future  Plans 

Numerical  analysis  of  normal  and  tangential  forebody 
blowdng  will  be  simulated  and  compared  to  the  asym¬ 
metric  strake  results  Additionally,  an  unstructured  grid 
approach  will  be  used  to  compare  Euler  and  Navier- 
Stokes  results  for  strake  and  blowing  applications. 

Publications 

Wurtzler.  K.  E.  "Numerical  Analysis  of  a  Chined  Fore¬ 
body  with  Asymmetric  Strakes."  AIAA  Paper  93-0051. 
January  1993. 
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Inviscid  Hypersonic  Flow  Around  a  Conically- Derived  Waverider 

Capt  Gregory  0-  Stecklein 
Interdisciplinary  and  Applied  CFD  Section 


Research  Objective 

The  current  high  level  of  interest  in  hypersonic  flight, 
brought  about  to  a  large  extent  by  programs  such  as  the 
National  Aerospace  Plane,  has  focused  much  attention 
on  the  stringent  requirements  of  aerodynamic  design  at 
these  demanding  flight  conditions.  The  major  objective 
of  this  research  is  to  analyze  the  inviscidly  dominated 
flowfield  of  a  waverider  forebody  for  on-  and  off-design 
flight  regimes. 

Approach 

The  governing  equations  are  the  three-dimensional  Eu¬ 
ler  equations  utilizing  a  finite-volume  solver  using  spa¬ 
tial  Roe-averaged,  flux-difference  splitting  They  are 
solved  with  a  .standard  McCormack  explicit  time  inte¬ 
gration  scheme.  Con\'ergence  is  accelerated  at  the  ex¬ 
pense  of  time  accuracy  by  the  local  time-stepping  tech¬ 
nique 

Accomplishments 

Three  general  classes  of  solutions  were  investigated  for 
two  levels  of  grid  refinement.  Inviscid  solutions  were 
obtained  for  on-design  conditions  of  Mach  10.  30.5  km 
pressure  altitude  and  zero  degrees  angle-of-attack  ( AoA  i 
and  two  off-design  flight  regimes.  Qualitative  results  il¬ 
lustrated  excellent  shock  capturing  and  flowfield  forma¬ 
tion  characteristics  with  grid  refinement. 


The  on-design  solutions  were  then  validated  against  the 
analytical  results  which  were  based  on  hypersonic  small 
disturbance  theory.  HSDT.  and  an  integral  momexitum 
development.  Quantitative  results  matched  well  with 
the  analytic  results  to  within  o'?!  for  lift  and  drag  co¬ 
efficients,  Off-design  calculations  were  run  for  the  wa¬ 
verider  at  Mach  10  for  an  Ao.A  range  of  ±5'T  and  at 
0“  AoA  and  a  Mach  range  of  6  to  15.  For  the  inviscid 
calculations,  the  off-design  perturbauons  were  charac¬ 
terized  by  wave  drag  effects  Of  particular  interest  was 
the  linear  wave  drag  bucket  associated  with  perturba¬ 
tions  to  the  on-design  attitude.  For  the  range  of  atti¬ 
tudes  mvestigated  in  this  research,  wave  drag  increased 
approximately  109c  per  degree  angle-of-attack. 

Future  Plans 

Current  plans  are  to  reevaluate  this  waverider  configu¬ 
ration  with  a  mass-addition  leading  edge.  The  implicit, 
full  N'avier-Stokes  mode  of  the  current  algorithm  will 
also  be  employed  to  investigate  the  effects  of  shock  wave 
and  boundary  layer  interaction 

Publications 

Stecklein,  G.  0  "Numerical  Solution  of  Inviscid  Hy¬ 
personic  Flow  Around  a  Conically-Derived  Waverider." 
AIAA  Paper  93-0320.  January  1993. 
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InviS'ud  Mach  Contours  .V/s-  =  10  Oio  =  0“ 


The  Flowfield  Past  the  X24C  Reentry  Vehicle 


Datta  Gaitond<='  -  I  niversal  Energy  Systems  Inc 
CFD  Researcli  Section 

Research  Objective 

The  purpose  of  this  research  is  to  examine  the  flowfield 
about  complete  lifting-body  configurations  at  hyper¬ 
sonic  speeds  utilizing  an  accurate  and  efficient  numeri¬ 
cal  algorithm.  The  ability  of  a  high  resolution  scheme  to 
provide  accurate  engineering  estimates  of  surface  ther¬ 
mal  and  mechanical  loading  is  examined. 

Approach 

The  full  3-D  compressible  Navier-Stokes  equations  in 
mass-averaged  variables  are  solved  in  a  finite  volume  for¬ 
mulation.  Roe's  flux-difference  split  method  is  used  for 
the  inyiscid  fluxes  while  viscous  fluxes  are  centered.  The 
implicit  numerical  algorithm  is  based  on  Gauss-Seidel 
line  relaxation.  Turbulence  closure  is  achieved  with  the 
Baldwin-Lomax  algebraic  eddy  viscosity  model. 

Accomplishments 

The  choice  of  algorithm  is  determined  with  a  crit¬ 
ical  examination  of  Ml'SCL  based  higher  order  up¬ 
wind  methods  and  a  range  of  limiters.  Spe¬ 
cific  emphasis  is  placed  on  the  prediction  of 
heat  transfer  rates  in  viscous  flows  The  flow 


past  the  X24C-iOD  configuration  is  then  solved  with  a 
mesh  of  roughly  half  million  grid  points.  The  so-called 
■'carbuncle'  anomaly  in  the  nose  region  is  eliminated  by 
a  combination  of  different  splittings  without  affecting 
the  overall  accuracy.  Analysis  of  the  computed  turbu¬ 
lent  flowfield  show  good  agreement  with  experimental 
surface  pressurf^  and  heat  transfer  data  Several  aspects 
of  the  computed  laminar  and  turbulent  flowfields  are  in¬ 
vestigated.  including  the  interaction  of  the  vortex  struc¬ 
ture  with  the  tail  fins 

Future  Plans 

The  theoretical  model' will  be  extended  to  incorporate 
high  temperature  effects  such  as  thermo-chemical  non- 
equilibrmn'i  and  ionization  A  ’2-0  version  of  the  code, 
employing  non-equihbnum  five-species  chemistry  and  a 
simple  harmonic  oscillator  model  for  vibrational  energy 
excitation,  has  already  been  developed 

Publications 

■Numerical  Simulation  of  Flow  Past  the  X24C  Reentr> 
X'ehicle,  '  D  Gaiionde  and  J  Shang.  AlAA  93-0319.  3l5t 
Aerospace  Sciences  Meeting  X  Exhibit.  Reno  1992. 
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Bomb  Blast  Simulation 

Wiliam  Z.  Strang 

Interdisciplinary  and  Applied  CFD  Section 
Research  Objective 

Simulate  air  blast  loads  induced  by  the  detonation  of  a 
small  charge. 

Approach 

The  general  approach  taken  weis  to  incorporate  a  bomb 
model  into  an  existing  computational  fluid  dynamics 
(CFD)  code.  COBALT  was  chosen  as  the  CFD  code  to 
host  the  bomb  model.  This  code,  an  unstructured  grid, 
higher-order  Riemann  solver,  was  chosen  becau.se  of  its 
robustness  and  the  ease  with  which  the  bomb  and  its 
surroundings  could  be  discretized.  The  bomb  model  is 
very  simple  owing  to  the  as.sumption  that  the  detailed 
physics  of  charge  detonation  are  unimportant  at  dis¬ 
tances  large  compared  to  the  charge  radius.  A  corollary 
to  this  assumption  is  that  the  mass  of  the  detonation 
by-products  is  small  with  respect  to  that  of  the  sur¬ 
rounding  medii  m.  With  these  assumptions,  a  charge 
may  be  simply  modelled  as  a  mass  and  energy  source. 
The  mass  and  energy  transfer  rate  is  assumed  constant 
and  is  calculated  from  tabulated  detonation  velocities 
for  the  charge  material  at  hand.  Owing  to  the  high 


inertial  forces  involved,  viscous  effects  were  not  simu¬ 
lated  Finally,  radiation  effects  were  assumed  negligible 
due  to  the  short  lime  scale  involved 

Accomplishments 

One  of  the  Aircraft  Loading  Tests  conducted  by 
WL/FI\'  was  simulated  In  the  test,  a  bare  spherical 
charge  was  located  at  the  centroid  of  a  rigid,  right- 
circular  cylinder  with  semi-oblate  spheroidal  endcaps. 
The  particular  test  simulated  was  chosen  for  its  accurate 
pressure  histories  at  several  recording  locations  The 
simulation  covered  the  first  11  milliseconds  of  the  blast 
event.  The  accompanying  figure  depicts  the  pressure 
contours  within  the  cylinder  1.56  milliseconds  after  det¬ 
onation.  In  it.  the  primary  blast  wave,  first  sidewall 
reflection  and  implosion  shock  are  seen  Comparisons 
between  experimental  and  numerical  pressure  histones 
at  various  recording  locations  show  excellent  agreement 

Future  Plans 

We  expect  to  simulate  a  detonation  tvuhin  a  rigid,  three- 
dimensional  fuselage 


Parallel  Processing  for  Computational  Fluid  Dynamics 

Stephen  Scherr 
CFD  Research  Section 


Research  Objective 

This  effort  describes  the  parallel  structure  of  an  explicit 
Navier-Stokes  algorithm  and  its  implementation  on  a 
massively  parallel  message-passing  computer. 
Approach 

The  analysis  of  complex  fluid  flows  and  the  future  cou¬ 
pling  of  CFD  techniques  in  interdisciplinary  applica¬ 
tions  require  ever  increeising  amounts  of  computational 
power.  Trends  in  computer  architecture  indicate  that 
the  required  computational  power  will  only  be  avail¬ 
able  through  the  use  of  many  processors  all  working  to¬ 
gether  to  perform  the  analysis.  The  architectures  which 
currently  show  the  most  promise  are  those  in  which 
each  processor  has  its  own  private  memory  and  com¬ 
municates  with  other  processors  by  psissing  information 
through  messages.  Samples  of  these  are  tightly  coupled 
massively  parallel  computers,  and  loosely  coupled  net¬ 
works  of  workstations. 

To  perform  a  computational  fluid  dynamics  analysis  in 
parallel,  the  entire  problem  must  be  divided  up  into 
pieces,  which  are  then  assigned  to  the  available  proces¬ 
sors.  For  best  performance,  this  decomposition  should 
evenly  distribute  the  work  among  the  processors,  and 
minimize  the  number  of  messages  which  need  to  be 
passed  between  processors.  In  addition,  support  mech¬ 
anisms  such  as  input/output  techniques  should  be  de¬ 
veloped  to  take  advantage  of  system  features 
Accomplishment  s 

An  existing  serial  implementation  of  the  two-step  ex¬ 
plicit  predictor-corrector  MacCormack  algorithm  for  so¬ 
lution  of  the  3-D  Navier-Stokes  equations  was  ported  to 


the  Intel  Touchstone  Delta  massively  parallel  computer 
The  parallel  design  features  a  2-D  pencilled  block  data 
decomposition  with  only  nearest  neighbor  communica¬ 
tions.  This  IS  well-suited  to  the  two-dimensional  mesh 
architecture  of  the  Delta  machine. 

As  a  test  of  the  implementation,  the  vortical  flow  around 
a  75®  swept  delta  wing  was  simulated  numerically  for 
flow  conditions  of  a  =;  30®.  Af-t  =  1.95.  and  a  Reynolds 
number  of  4.48  x  10®  .  The  computational  grid  was  an 
80  X  80  X  80  mesh  with  an  H-H  topology.  Timing  re¬ 
sults  for  the  main  computational  loop  indicate  a  factor 
of  10  speedup  over  a  single  Cray-2  proces.sor  when  oOO-i- 
processors  of  the  Delta  machine  are  utilized  (see  chart 
below).  The  best  performance  occurs  when  the  work 
load  is  evenly  balanced  across  the  processors.  Various 
methods  for  remapping  a  13  x  39  set  of  blocks  into  a 
16  X  32  physical  processor  layout  were  examined.  The 
best  layout  occurs  by  maximizing  the  number  of  pro¬ 
cessors  having  nearest  neighbor  communications.  This 
configuration  is  illustrated  in  the  diagram  below. 
Future  Plans 

Future  work  in  parade'  rocessing  for  computational 
fluid  dynamics  will  include  implementation  of  a  stan¬ 
dard  turbulence  model  under  a  blocked  data  decompo¬ 
sition.  extension  of  techniques  to  other  algorithms,  such 
as  an  explicit-implicit  flux-split  algorithm,  and  investi¬ 
gation  of  methods  for  parallelizing  unstructured  grids. 
Publications 

Scherr.  S.3..  “Implementation  of  an  Explicit  Navier- 
Stokes  Algorithm  on  a  Distributed  Memory  Parallel 
Computer."  AlAA  Paper  93-0063.  1993 
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Advanced  Technology  Air  Vehicle  Geometry  Data  Base 

Howard  T  Emsley 

Interdisciplinary  and  Applied  CFD  Section 


Research  Objective 

Accurate  geometric  representations  of  USAF  aircraft 
and  weapons  are  necessary  to  perform  many  Air  Force 
support  activities  such  as  weapon  and  vehicle  modifi¬ 
cations.  Computational  Fluid  Dynamic  (CFD)  analysis, 
and  weapon  fit  checks.  Acquiring  necessary  geometries 
is  a  difficult  task  due  to:  a  lack  of  geometric  standard¬ 
ization  in  the  industry,  the  proliferation  of  computer- 
aided  design  tools,  and  the  absence  of  complete  L’S.AF 
aircraft  external  geometries.  Standardization  of  geome¬ 
try  formats  under  the  Initial  Graphics  Exchange  Spec¬ 
ification  (IGES).  Product  Data  Exchange  using  STEP 
fPDES).  and  Computer-aided  Acquisition  and  Logistic 
Support  (CALS)  jTiay  be  the  first  step  toward  resolv¬ 
ing  the  problem,  but  another  issue  involves  the  acquisi¬ 
tion  and  storage  of  air  vehicle  geometries  for  use  by  Air 
Force  organizations.  With  accurate  geometries  in  hand 
the  timeliness  and  ability  to  perform  support  activities 
can  be  improved,  duplication  of  effort  can  be  reduced, 
and  common  geometric  configurations  can  be  provided 
to  Air  Force  research  and  support  organizations. 

Approach 

To  perform  support  activities  for  Air  Force  systems,  ge¬ 
ometries  of  USAF  aircraft  and  vveapons  should  be  ac¬ 
quired  as  advanced  surface  representations  and  stored 
in  a  data  base.  The  goal  of  this  program  is  to  reduce 
the  time  required  to  perform  CFD  analyses  by  providing 
accurate  geometric  representations  to  CFD  users. 


To  reach  our  goal,  a  data  base  structure,  a  set  of  geome¬ 
try  formats,  and  the  tools  to  store,  retrieve  and  interac¬ 
tively  view  complete  air  vehicle  geometries  are  required. 
Concurrent  with  the  development  of  the  system,  WL  is 
collecting  current  and  future  USAF  air  vehicle  geome- 
trie.s  to  store  in  the  data  base.  Data  formats  specified 
by  Cals  will  be  used  to  acquire  the  geometries,  and 
MIL  standard  modifications  or  clarifications  will  be  ad¬ 
dressed  to  reduce  problems  in  translating  geometric  rep¬ 
resentations  from  different  systems 

Accomplishments 

-A  contract  to  develop  the  data  base  structure  and  inter¬ 
face  software  was  issued  to  Rensselaer  Design  Research 
Center  in  May  1992.  A  non-uniform  rational  b-spline 
geometry  format  (STEP  NURBS)  was  selected  and  the 
preliminary  interface  software  has  been  developed  Ad¬ 
ditional  work  in  progress  involves  interfacing  with  ex¬ 
isting  standards,  and  the  inclusion  of  topology  between 
geometric  entities.  The  former  involves  format  conver¬ 
sion  of  IGES  entities  into  STEP  NURBS  entities  for 
storage  and  retrieval. 

Future  Plans 

As  the  contracted  effort  continues,  air  vehicle  geometries 
will  be  collected  for  conversion  and  storage  in  the  data 
base.  In-house  efforts  associated  with  this  effort  are  con¬ 
centrating  on  efficient  ways  of  interfacing  .NURBS  sur¬ 
faces  with  existing  and  developing  grid  generation  tools 
for  CFD. 


.'\d vatii'e<l  Technology  .\ir  Vehicle  G'>ometry  D.it.i  Base  Browser  and  \  lewer 
fl’onzontal  tail  data  i.s  displayed  in  brow.=;er  .and  highlighted  in  viewer) 
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Computational  Electromagnetics 

J.  S.  Shang 

Flight  Dynamics  Directorate 
Research  Objective 

The  ultimate  objective  of  the  present  effort  attempts 
to  establish  a  numerical  simulation  capability  of  elec¬ 
tromagnetic  phenomena  for  air  vehicle  technology.  The 
basic  process  involves  the  transition  of  innovative  tech¬ 
niques  from  computational  fluid  dynamics  (CFD)  to 
computational  electromagnetics  (CEM)  for  signature 
analysis  and  control,  target  identification,  microwave 
communication,  as  well  as  environmental  hazard  pro¬ 
tection.  Specific  objectives  are  to  first  develop  efficient 
numerical  procedures  for  solving  time-domain,  three- 
dimensional  Maxwell  equations,  and  then  map  numeri¬ 
cal  algorithms  onto  massively  parallel  processors. 

Approach 

The  time-dependent  Maxwell  equations  for  propagating, 
refracting  and  diffracting  eletromagnetic  waves  consti¬ 
tutes  an  initial  value  problem.  The  hyperbolic  system 
is  solved  by  characteristic-based  schemes,  which  are  de¬ 
rived  from  eigenvalue  and  eigenvector  analysis.  The 
three-dimensional  system  of  equations,  when  written  in 
flux  vector  form,  can  be  di  .gonalized  in  each  single  di¬ 
mension,  The  concept  of  flux-difference  splitting  will 
permit  windward  schemes  to  mimic  the  wave  mechanism 
of  information  propagation.  In  three-dimensional  appli¬ 
cations.  the  cyclic  fractional-step  or  the  time-splitting 
formulation  ensures  a  .second-order  accurate  procedure 
in  time  and  space.  Both  finite-difference  and  finite-  vol¬ 
ume  solving  procedures  will  be  implemented  on  a  gen¬ 
eral  curvilinear  coordinate  system  to  accommodate  com¬ 
plex  scatter  shapes.  Extremely  high  numerical  efficiency 
for  high  frequency  wave  simulation  are  being  achieved 
through  a  domain  decomposition  technique  to  minimize 
the  data  communication  on  a  massively  parallel  proces¬ 
sor. 


Accomplishments 

The  characteristic-based  methods  have  been  validated 
by  applying  them  to  a  three-dimensional  wave  guide  and 
a  singular  electric  dipole  problem  The  present  approach 
has  demonstrated  a  unique  potential  in  alleviating,  and 
under  some  special  circumstances  m  eliminating,  the 
fundamental  dilemma  of  reflecting  waves  from  the  trun¬ 
cated  spatial  domain.  The  numerical  procedures  also 
exhibit  significantly  improved  numerical  .stability  prop¬ 
erties  and  reduced  alaising  errors  from  the  current  com¬ 
putational  electromagnetic  methods  in  use  .4  very  high 
data  processing  rate  of  six  gigaflops  has  also  been  at¬ 
tained  by  a  present  finite-difference  procedure.  Most 
importantly,  on  a  massively  parallel  system,  the  present 
second-order  accurate  finite-difference  procedure  is  scal¬ 
able.  The  scalability  is  sustained  up  to  512  nodes,  with 
only  a  3  percent  parallel  efficiency  degradation. 

Future  Plans 

Further  research  will  be  pursued  to  refine  the  finite- 
volume  procedure  for  an  air  vehicle  radar  signature  anal¬ 
ysis  effort.  A  comparative  study  is  also  planned  to  as¬ 
sess  the  relative  merits  of  the  finite-volume  and  finite- 
difference  approaches  in  application  to  complex  scatter 
geometries,  and  numerical  efficiency  on  massively  par¬ 
allel  computers.  The  characteristic  concept  will  finally 
be  extended  to  develop  a  meshless  computing  procedure 
for  efficient  and  accurate  CEM  applications  and  diverse 
dual-use  purposes 

Publications 

Shang.  J.  S..  '‘A  Fractional-Step  Method  for  Solving  3- 
D  Time  Domain  Maxwell  Equations."  AiAA  paper  93- 
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